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Abstract. We are conducting a project aimed at surveys and repeated observations of 
red variables (or long-period variables) in globular clusters. Using the IRSF/SIRIUS near- 
infrared facility located at South Africa, we are observing 145 globular clusters that are 
accessible from the site. In this contribution, we present our observations and preliminary 
results. We have discovered many red variables, especially in the Bulge region, whose mem- 
berships to the clusters remain to be confirmed. Using a sample of all red variables (both 
already known and newly discovered ones) in globular clusters except those projected to 
the Bulge region, we produce a log P-K diagram and compare it with those for the Bulge 
and the Large Magellanic Cloud. A prominent feature is that the bright part of overtone- 
pulsators' sequence (B + and C) is absent. We discuss its implication on the evolution of red 
variables. 
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1. Introduction 

The study o n red variables ha s advanced very 
much since IWood et al. I (Il998l) found paral- 
lel sequences on the period-magnitude dia- 
gram as a result of the MACHO project for 
the Large Magellanic Cloud. Many papers 
have been published on the variables in the 
Magellanic Clouds (Ita et al. 2004a, 2004b; 
Kiss & Bedding, 2003, 2004; also see the con- 
tribution by Dr. Kiss in this volume). However, 
it is difficult to discuss the stellar evolution di- 
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rectly because the Magellanic Clouds contain 
red variables with different ages and different 
chemical components. 

In contrast, globular clusters (GCs) are 
systems of single stellar population, so that 
we can tell the age and the initial chemi- 
cal component for each variable star. Making 
use of this advantage, some authors have 
worke d on the e volution of red variables in 
GCs. [yVhitelock] jl986l) found Mira variables 
and semi-regular variables in GCs define a 
relatively tight sequence on the log P-Mboi 
diagram. iBedding & Ziilstra I dl998l) consid- 
ered the sequence as an evolutionary track, 
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whose slope is c onsistent with an evo lution- 
ary calculation by Vassiliadis & Wood ( 1993). 
On the other hand. iPerl & Tuchmanl ( 1 19901) 
theoretically examined the logP-Mboi dia- 
gram with a staircase pattern for variables 
in 47 Tuc, and concluded that a switch of 
pulsational- mode occurs at a certai n luminos- 
ity. Recently, Lebze her et al. I d2005l) also pro- 
posed the switch of mode based on log P-K di- 
agram and their measurements of variations of 
radial velocities. 

These works clearly demonstrated that 
GCs are good testbeds for studying the evolu- 
tion of red variables. However, early surveys 
were not complete, and many of them were 
optimized to search for short- period var iables 
such as RR Lyrae variables dClement et al. I 
1200 ll) . In addition, they are mostly conducted 
in visible wavelength. Some red variables 
are undergoing heavy mass loss and holding 
thick circumstellar materials, w hich make the 
stars only detectable in infrared! Tanabe et al. 
Il997l) . So that we decided to conduct a near- 
infrared survey for red variables in GCs. 

2. Observations and Preliminary 
Results 

We started our near-infrared survey in the 
spring of 2002, using the IRSF 1.4 m 
telescope and the SIRIUS camera. They 
are established by Nagoya University and 
National Astronomical Observatory of Japan, 
and sited at Sutherland station of South African 
Astronomical Observatory. They can take three 
images in J-, H-, and /T s -band simultane- 
ously with the 8' x 8' fields-of-view. For 
the details of the IRSF and the SIRIUS, see 
Naeashima et al. I (Il999l) andtNagavam a et aT71 

There are 150 GCs f ound in the Mil ky Way 
and the neighborhood (Harris 1996). From 
the IRSF site, we can reach to about +30° in 
Declination, and we have observed such 145 
clusters at least once. Until the end of 2004, we 
have observed 133 clusters more than 10 times 
and 40 clusters more than 20 times. 

With data analyses for the 133 clusters, we 
detected variations of 58 known red variables 
in 30 clusters and discovered 185 red vari- 



ables in 44 clusters. Detection limits of vari- 
ability is typically better than 0.1 mag. Many 
of the newly discovered variables lie in the di- 
rection of the Galactic Bulge. 31 clusters with 
the new variables locate within the \l\ < 10° 
and \b\ < 10° region. Because the Bulge itself 
contains many red variables, it is important to 
check their memberships. 

Especially, it is urgent for some peculiar 
objects. Among the new ones, there are very 
red stars, e.g. (J - K)q > 3. It is not known 
that such stars exist in GCs. We detected SiO 
maser emissions, evidences of relatively thick 
circumstellar materials, from four of the newly 
discovere d variables and also from two known 
variables (Matsunag a et al. 1200 5a). The veloc- 
ities of the emissions indicate five of them can 
be associated with GCs. 

Using the near-infrared surface bright- 
ness obtained by the COBE/DIRBE data 
( Matsunag a~et al l2005bll . roughly 20 % of the 
detected variables are considered to be associ- 
ated with the GCs. 

3. Discussions 

3.1. The log P-K diagram and 
comparison with the LMC 

In the following, we will discuss the log P-K 
diagram using our samples. In order to avoid 
the Bulge field contamination, we adopt vari- 
ables whose galactic coordinates are out of 
the Bulge range: |/| > 10° and/or \b\ > 10°. 
Collecting red variables with known periods 
and near-infrared measurements, our sample 
contains 86 red variables in 29 GCs. For 27 
stars among them, we detected their variations 
and obtained their mean magnitudes (mean of 
minimum and maximum); we did not detected 
significant variations for 50 stars; 9 stars are 
located too north to be observed and we just 
adopted the magnitudes listed in the 2MASS 
point source catalog. 

Fig. [l] shows the log P-K diagram. 
Filled circles for variables in metal-rich 
GCs ([Fe/H]> -1) and triangles for those 
in metal-poor ones are superposed on gray 
points for those in th e Large Magellanic 
Cloud dlta et al. I l2004bl) . Mira variables in 
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Fig. 1. The log P-K diagram. Filled circles for variables in metal-rich clusters ([Fe/H] > -1), star symbol s 
for those in metal-poor ones, and gray points for the data of the Large Magellanic Cloud ftta et al. l2004bl) . 



the metal-rich clusters obey the sequence C, 
and some variables are located at around the 
sequence B (notati ons for the se quences are 
based on those of lltaetal.l l2004a). As is well 
known, bright Mira variables belong to o nly 
metal-rich GCs dFroeel & Whitelockll99l . 

Few variables in GCs lie on the sequences 
A, B + , C", and D. For variables on the se- 
quences A and D, the periods are difficult to 
get. Those on the sequence A have small am- 
plitude and irregular light curves, and ones 
on the sequence D have very long periods 
and sometim es double-periodic light curves 
( Wood1 l2000l) . We need massive photometric 
data and very careful analyses in order to con- 
clude the absence of these sequences. 

On the other hand, variables on the se- 
quences B + and C should be easier to find, 
at least compared with those on the sequence 
B . It is also known about GCs that all the 
stars brighter than the tip of Red Giant Branch 
(RGB ) are found to be Mira variables (see 
iFrogel & Whitelock 1119981 for example). So 
that, we conclude this absence is a real fea- 
ture of variables in GCs. Please note that this 



absence was ut ilized to propose the evolu- 
tionary track by iPerl & Tuchman I ill 9901) and 
Lebzelteret al. (2005). We confirmed the ab- 
sence from a larger sample including many 
GCs. 



3.2. Comparison with the Bulge and its 
implication 

We found the intermediate part of the sequence 
C'(or B + ) exists on the log P-K diagram for the 
Galactic Bulge given by Glass & Schultheis 
(2003; see their Fig. 7). On that sequence, there 
is an vacant region at the brightest part, but 
some variables are as bright as Mk = -7.5 with 
the distance modulus assumed to be 14.5 mag. 
Although it is difficult to distinguish between 
the sequences C and B + , they are bright over- 
tone pulsators that are not found in GCs. The 
top of the sequence is brighter than the tip of 
RGB. Although there are some uncertainties 
about the tip of RGB for the Bulge, such as 
spatial depth and metallicity effects, we con- 
sider they should work on both the sequences 
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A and B. We conclude the intermediate part is 
in the Asymptotic Giant Branch (AGB) phase. 

It is suggested that AGB stars in the Bulge 
take the evolutionary track with the staircase 
pattern similar to that of 47 Tuc, but the 
luminosity at which the mode of pulsation 
switches is different. This can be interpreted 
that the variables in the Bulge have larger ini- 
tial masses than those in 47 Tuc and take 
a brighter evolutionary trac k as presented by 
Vassiliadis & Wood I lll993l) . The star forma- 
tion history of the Bulge is still controversial, 
but the majority of the constitu ent stars have 
the si milar ages to those of GCs llZoccali et alJ 
120031) . Nonetheless, their high metallicities 
have an effect on initial masses of the red vari- 
ables. Theoretically, red giants with the high 
initial metallicities are expected to have the 
larger initial masses (Iben & Renzini 1984). 

It is important to discuss GCs in the Bulge 
region whose metallicities are similar to the 
Bulge field population after confirming the 
memberships of variable stars. It is also inter- 
esting to investigate the log P-K diagram for 
some massive clusters with the intermediate- 
age (~1 Gyr, i.e. more massive) AGB stars in 
the Magellanic Clouds. 

4. Conclusions 

We discovered many red variables in the course 
of our near-infrared survey toward 145 GCs. 
Many of the newly discovered variables lie in 
the direction of the Galactic Bulge. Since there 
is a large field population of red variables in the 
Bulge, it is important to check their member- 
ships by radial velocities and proper motions. 

Using the sample outside the Bulge region 
including both the new variables and already- 
known ones (86 variables in 29 clusters), we 
discussed the distribution on the log P-K dia- 
gram. We found the absence of the bright part 
of the sequences C and B + , which was known 
for 47 Tuc, with more variables whose metal- 
licities are widely spread. We also compared 
the distribution with that of the Bulge. It is sug- 
gested that AGB variables switch from over- 
tone pulsators to Mira variables at a certain lu- 
minosity, which depends on stellar parameters 
such as initial mass and metallicity. 
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